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Abstract 

Bacteria are able to respond to environmental signals by changing their rules of movement. 
When we take into account chemical signals in the environment, this behaviour is often 
called chemotaxis. At the individual-level, chcmotaxis consists of several steps. First, the 
cell detects the extracellular signal using receptors on its membrane. Then, the cell processes 
the signal information through the intracellular signal transduction network, and finally it 
responds by altering its motile behaviour accordingly. At the population level, chemotaxis can 
lead to aggregation of bacteria, travelling waves or pattern formation, and the important task 
is to explain the population- level behaviour in terms of individual-based models. It has been 
previously shown that the transport equation framework |12[ 113) is suitable for connecting 
different levels of modelling of bacterial chemotaxis. In this paper, we couple the transport 
equation for bacteria with the (parabolic/elliptic) equation for the extracellular signals. We 
prove global existence of solutions for the general hyperbolic chemotaxis models of cells which 
process the information about the extracellular signal through the intracellular biochemical 
network and interact by altering the extracellular signal as well. The conditions for global 
existence in terms of the properties of the signal transduction model are given. 

1 Introduction 

The flagellated bacteria (e.g. Escherichia coli, Salmonella typhimurium, Bacillus subtilis) are 
single-celled organisms. They are usually too small to be visible by the naked eye; typically, 
they have the size of microns (see |25| ll8j for review). The behaviour of a bacterium is primarily 
influenced by concentrations of various chemicals inside the cell. Since bacteria are small, we can 
assume that the concentrations of the chemicals inside the cytoplasm are uniform. Therefore, 
we can suppose that the cells are points. Moreover, to create a mathematical description of a 
bacterium, we introduce the vector of internal state variables jl21 1131 llUj 

y = {yuy2,...,ymf em^, (1) 

where yi, i = l,...,m, are concentrations of various chemicals (proteins, receptor states etc.) 
inside the cell involved in the processes of interest. The individual behaviour of a cell primarily 
depends on the vector y which is a function of time. Consequently, the state of bacterium is 
uniquely determined by the vector (i, x, v, y) where x G is the position of a cell, v G is 
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its velocity, y E M'" is its internal state, t is time and = 1,2, or 3, is the dimension of the 
physical space. 

During its life, a cell must communicate with its environment in order to find nutrients, to 
avoid repellents, to find mates etc. For this purpose, there are receptors in the cellular membrane 
which can detect various chemicals in the environment. We describe the chemicals outside the 
bacterium by the signaling vector (which depends on the position of the cell x and time t) 

S(x,t) = (5i,52,...,5A/f GM^^ (2) 

Then the evolution of the internal state vector y depends also on the signaling vector S. Since 
we describe chemical processes, we can assume that y evolves according to system of ordinary 
differential equations 

^=F(S(x),y). (3) 

This system formally captures all biochemistry inside the cell and therefore, the concrete form 
of the vector function F : M*^ x M™ can be very complicated depending on the number of 

details which are included in the model. 

Bacterial movement and the signal transduction network @ will be discussed in more details 
in Section |3J From the mathematical point of view, the movement of the flagellated bacteria can 
be viewed as a biased random walk. The properties of this random walk depend on the internal 
state y and bacterial velocity v. The classical description of the bacterial movement is the so 
called velocity jump process [1011^ El- It means that the bacterium runs with some velocity 
and at random instants of time it changes its velocity according to the Poisson process with the 
intensity A(y). 

Let /(x, v,y,t) be the density function of bacteria in a (2A^ + m)— dimensional phase space 
with coordinates (x,v,y) where x G is the position of a cell, v G F C is its velocity 
and y G is its internal state, which evolves according to ©. Thus /(x, v, y, t)dxdvdy is the 
number of cells with position between x and x + dx, velocity between v and v + dv, and internal 
state between y and y + dy. Then evolution of / is governed by the following transport equation 



dt 



+ Vx • v/ + Vy • F(S(x), y)/ = -A(y)/ + / A(y)i^(v, v', y)/(x, v', y, t)dv' (4) 

Jv 



where the kernel ^^(v, v',y) gives the probability of a change in velocity from v' to v, given 
that a reorientation occurs. We assume that the random velocity changes are the result of a 
Poisson process of intensity A(y) 0. The kernel K is non-negative and satisfies the normalization 
condition 

/ i^(v,v',y)dv = l (5) 
Jv 

where 

y is a symmetric compact set in M^. (6) 

Realistic examples of the kernel K{v, v', y), set V, signal transduction network F and the turning 
frequency A(y) are given in Section |3I They all satisfy the above basic assumptions. 

To write equation ^ in more compact form, we introduce the kernel T defined as a product 
of the turning frequency A and the kernel K, i.e. 

r(v,v',y) = A(y)K(v,v',y). (7) 
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Moreover, our goal is to couple equation (jlj with the realistic system of partial differential 
equations for the extracellular signal vector S. We assume that the external signal diffuses. It 
can be also produced or degraded by bacteria, degraded on its own or the components of S can 
react with each other in the extracellular space. Hence, the general hyperbolic system of interest 
can be written in the following form: 



where n = n(x, t) is the macroscopic density of individuals at point x G and time t given as 



D is a diagonal M x M matrix which diagonal elements are diffusion constants of different 
chemicals in the extracellular signal vector S and the term R : x R ^ describes the 
creation, reaction and degradation of the signals. 

The goal of this paper is to prove global existence results for the system (jSJ - ©. We will 
focus on one-dimensional case in what follows. In Section [21 we will start with a simple model of 
the signal transduction (|TT|) which was used previously ^L2i tl^ • The simple model (fTTj) has the 
essential properties of the realistic models of the signal transduction, but it is more tractable from 
the mathematical point of view than more complex models of bacterial chemotaxis. We prove 
the global existence of solutions of the one-dimensional version of - © with the simplified 
model of the signal transduction. 

In order to study the general case, we first review the relevant biology in Section |31 This will 
help us to specify the realistic conditions on signal transduction model turning frequency 
A(y), turning kernel -ftr(v,v',y), set V, diffusion matrix D and the reaction term R(S,n) in 
equations (jSJ -© • In Section ^ we study the global existence for the system (jSJ - ® for 
general models of bacterial signal transduction which are introduced in Section |21 We also 
consider that equation © is at quasi-equilibrium, i.e. we consider the elliptic equation for the 
signal in Sectional 

Hence, this paper consists of two main mathematical results. First, we prove the global 
existence of solutions to the problem (jH)) - (jU for the simplified model of signal transduction 
(|TT|) and for the system of parabolic equations Q for the extracellular signal (see Section |21 
Theorem^. Then, we prove the global existence of solutions for the general model of signal 
transduction © coupled with the system of elliptic equations for the extracellular signal (see 
Section!^ Theorem The necessary growth assumptions on turning frequency T(v,v',y) 
are given in terms of the signal derivative along the cell trajectory. It means that the growth 
estimates on T include the temporal derivative as well as the spatial derivative of the extracellular 
signal. Finally, we provide discussion and comparison with relevant results from the literature 
in Section |S1 




dS 

at 



L>AS + R(S,n) 



(9) 




(10) 
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2 Global existence for a simplified model of signal transduction 



A simplified model of excitation-adaptation dynamics was studied in ^1 ^1 |^ where y 
(2/1)2/2)"^ S I^^ and the right hand side of equation (jHJ was given as 



F2 



/ 5(S(x,t))-(yi + y2) \ 

V 



c/(S(x,t)) - y2 



(11) 



ta 



where and ta are positive constants and g : [0,oo). We will see in Section that 

the simplified model has the essential properties of realistic signal transduction models. 
Hence, the model pij) is a natural starting point of this paper. For simplicity, we work in a 
one-dimensional physical space, i.e. = 1, and the goal of this section is to prove Theorem ^ 
about the system (jSJ ~ ©• In what follows, we denote L^{^), 1 < p < oo, Q C M'^, the Banach 
space of measurable functions with the finite norms 



I LP 



|/i(x)|Pdx 



1/p 



for 1 < p < oo. 



and 



ess sup \ h(x) 
n 



We denote W''^P{n), 1 < p < oo, n C R'^, the usual Sobolev space 
W'^'PiCt) = \h£ LP{Vt) I Mcx £ \a\ < k ^'"'^ 



...x„ 



where ot = (01,021 ■ ■ ■ , ctd) G Ng is a vector of nonnegative integers and \a\ = oi + 02 + • • • + ad- 
The norm in is defined as 



KyJu-^ 2 * * * d 



LP 



To simplify mathematical formulas, we will make use of the following notation. For any function 
/i : M'^ ^ M, V/i denotes the gradient of h with respect to all variables and Vx^x2h is the 
2-dimensional gradient vector with respect to the variables xi and X2 only, i.e. 



V/i 



dh dh dh \ 
dxi ' dx2 ' ' dxd J 



and 



dh dh 

dxi ' dx2 



(12) 



We already made use of this notation in equation Q where the gradients of the function / were 
taken only with respect to the selected parts of the state vector. In this section, we study the 
movement of cells in one dimension, i.e. N = 1. Moreover, we assume that the external signal 
diffuses and it is produced by bacteria and degraded on its own. Hence, the system of equations 
(jEJ ~ © reads as follows 



dt 



+ V.-^/ + VyF(S(x),y)/ 



V 



~dt ~ ~dx^ 



T{v,v',y) f{v') - f{v) 
i = l,...,M, 



dv' 



(13) 
(14) 
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where di, ki and are positive constants and n = n{x, t) is the macroscopic density of individuals 
at point x G M and time t given by (|1U() . Position x and velocity v are scalars for = 1, so we 
do not use bold letters for position and velocity in equation ()13() . Otherwise, equation (|13() is 
the same as equation (jSJ. Following notation (|12jl . symbol V^/ denotes the partial derivative of 
distribution function / with respect to x. Let us note that (depending on the form of function g 
in some extracellular signals might be attractants and some extracellular signals might be 
repellents. If we have sufficient growth estimates on function g and kernel T, we can guarantee 
the global existence of solutions of system - (|14|) as it is shown in the following theorem. 

Theorem 1 Consider that the function F is given by (| Assume that there exist non- 
decreasing positive continuous functions £ C (R) satisfying 



|g(z)| + |V<7(z)| <$(|z|) and \T{v,v' ,y)\ + \VT{v,v' ,y)\ < ^ {\y\) . 



(15) 



Assume that /o G VF^'^(M x V x M?) n W'^'°°{R x V x M^) with compact support and So G 
[147i'°°(M)]*^ with compact support. Then there exist global solutions of the system (ji^p - ^141 
satisfying 

f{; •, •, t) e W^'^{R X y X M^) n W^'°^{R xV x M^), (16) 



and initial conditions /(•, •, •, 0) = /o(-, •, •) and S(-, 0) = So(-). 



(17) 



First, the characteristics of the hyperbolic equation H13() are given for = 1 as 

^ = V, ^ = 0, ^ = F(S(X(s),.),Y(.)). 

as as as 

Then along back-time characteristics starting at {x,v,y,t), we have for < s <t, 
X {s]x,v,y,t) = X — V {t — s) , 
Y{s;x,v,y,t) = y- T F (S (X(t), r) , Y(t)) dr. 

J s 

Next, we will prove several auxiliary lemmas. 



(18) 

(19) 
(20) 



Lemma 1 Derivation of the characteristics (| j,9|) and ^20fi with respect to the initial conditions 
gives, for < s <t, 



dX dY 
ox ay 



dF 

dy 



(s-t) 



where 



dF 

dy 



1 

~t~a ) 







(21) 



Moreover, 



, dY 

det — — = exp 
oy 



^ + f )(^--) 



> 1. 



(22) 
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Proof. We differentiate (PU)) with respect to y to get 

^ - I + r dF dY 
dy ^ Jt dy dy 

wliere I2 is tlie 2x2 identity matrix. Let 
then we have 



Using we obtain 



(9F (9F 



Integrating the last equation, we have 

G (s) = exp 



dF 
dy 



(s-t) 



(23) 



which deduce (j2T|). Computing the determinant of (|2T|1 . we derive (j22|). 



Lemma 2 Lei us assume anc? (| T/ien i/ie solution of (\18^ satisfies 



Q.E.D. 



|Y(r)| <C{l + ^{ sup |S(X(s),s) 

I VO<S<T 

where C depends on the y-support of fo and So, ta, and t^. 



(24) 



Proof. Using the assumption (|15)) and applying the Gronwall inequality to the ordinary differ- 
ential equation ((TH)) yields 

Y2{t) = Y2{0)exp{-T/ta) + ^ r g{S{X{s),s))eMis-T)/ta]ds 

ta Jo 



< |y2(0)| + ^$( sup \S{X{s),s) 

*a V0<s<r 



In a similar way, we get 



1 



yi(r) = n (0) exp (-r/te) + - / {g {S^X (s), s)) - {s)} exp[is - T)/te]ds 

te Jo 

< |yi(0)| + i sup \Y2is)\ + ^<^( sup \SiX{s),s) 

te 0<s<T te \0<s<T 



Thus we deduce (^1)) . 



Q.E.D. 



Lemma 3 If n £ L°°{[0,oo) : (M) H (IR)), then the solution S of the system of equations 
di^D satisfies 



as , , 



l|S 

< c 



L^<C sup ||n(t)||^i =C||n (0)11^1, 

0<T<t 



ln( sup ||n(r)||^2) 

VO<T<t / 



l + ll^(0)ILi (^l + (lnt)+ + 
where (•)_!_ means the positive part and the constant C depends only on ki, and d. 
Proof. See ^1 Lemma 4]. 



Q.E.D. 



Proof of Theorem ^ Integrating ((Tl^ along the characteristic ((T^ - ^Z0\ from to t and 
using (|T5|l . we get 

/(x,T;,y,i) </o(X(0),^,Y(0)) + 

ft 



+ 



Cf *(|Y(t)|)x f{X{T),V,Y{T),T)+ [ f{X{T),v',Y{T),T)dv' 

Jo I Jv 

+ r I Vy • F (S(X (r)), Y (t))| / (X (r) , v, Y (r) , r) dr. 
Jo 



dT+ 



Since • F = — ^ — ^, we get (using Lemma |2)) 



f{x,v,y,t) < fo{X{0),v,Y{0)) + C [' fiX{T),v,Y{T),T)dT+ 

Jo 



(25) 



+C ^ C 



! + $( sup |S(X(s),s)| 

.0<s<r 



\V\fiXir),v,YiT),T)+ [ f{X{T),v',YiT),r)dv' 

Jv 



dr 



where C is a constant depending only on support of /o, 5*0 , te and ta- Using Lemma Q we have 



dY\ 



-1 



(26) 



and so 



[ [ fP{X{T),v',Y{T),T)dv'dxdvdy 

JRxVxRM Jy 



= \V\JfP{X{T),v',Y{T),T)(^det^^ ' [det^^ 'dv'dXdY < 

< \V\ J fP {X (r) , v', Y (r) , r) d^^'dXdY. 
Taking the p-th power of ()25() and integrating over x, v, and y yields 



(27) 
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l + $( sup S(X(s),s) 

.0<s<t 



X / ||/(r)||^,dT. 



Lemma |21 implies 



sup ||S(-,r)||^^ <C sup ||n(t)||^, =C7||n(0)||^, <C||/o||^,. 

0<T<t 0<T<i 

Consequently, using (|27|) and (|28|) . we obtain 

ll/WII^P < ll/oll^. + C{l + ^{C[l + $(C||/o||^J])} X 1^ ||/(r)||^,dr. 

Applying the Gronwall inequality, we obtain, for all 1 < p < oo, 

< C(fci,A:0,di,te,ta,||/o||^p,supp/o,^',$,|y|) <oo. 



(28) 



(29) 



(30) 



We now compute a priori estimates on derivatives of /. We differentiate H13|) with respect of x, 
integrate along the characteristic ((T^ - (PU]) from to t and use ((T3|l to get 



^(x,T.,y,i) 



< 



||(X(0),«,Y(0)) 



+ C / ^(|Y(r)|)x 



1^1 



|^(X(r),.,Y(r),r) 



+ 



|^(X(r),.',Y(r),r) 



df' 



dT+ 



+ /VyF(S(X(r)),Y(T)) 

JO 



+ 



'0 



5/ 
9x 



(X(r),r;,Y(T),r) 



dT+ 







|Vy/(X(T),»,Y(T),T)|dT. 



Similarly, differentiating ()13() with respect of yi or 2/2, integrating along the characteristic (|19|1 - 
(Pn|) from to t and using ((T^. we obtain 

|Vy/(x,i;,y,t)| < |Vy/o(X(0),^;,Y(0))|+C r^(|Y(X(T),r)|)x 

^0 

X { 1^1 [ I/I + I Vy/I ] (X(r), t;, Y(r), r) +l^[\f\ + I Vy/| ] (X(r), i;', Y(r), r) dv'}dT+ 

+C r I Vy • F (S(X (r)), Y (r)) 1 1 Vy/ (X (r) , Y (r) , r) |dr. 
Jo 

If the interior of set V is nonempty, we can also define the derivatives of / with respect of v for 
any point in the interior of set V. Differentiating (|13() with respect of v and integrating along the 
characteristic ()19p - ()2U() from to t, it implies 



■^{x,v,y,t) 



< 



^{X{0),v,Y{0)) 



+ 







dr+ 



+ C ^(|Y(r)|)x |y||/(X(T),t;,Y(r),r)| + 
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+1^1 



dl 

dv 
t 



+ jjf{X{T),v',Y{T),T)\dv'Yr+ 



{X{t),v,Y{t),t) 

+ f\VyF{S{X{r)),Y{r)) 
Jo 

Using (jHOJ, (P^. Lemma 121 and Gronwall inequality, we deduce 



^{Xir),v,Y{r),r] 



dr. 



9// N 



+ 



LP 



+ ||Vy/(t)||^, < 



(31) 



LP 



< C [ki,k^,di,te,ta, ll/oll^^ip, ||S (0) 11^^ j,,supp/o,suppSo,^,^>, \V\j < oo. 
Combining (|30|) and pT|) . we obtain (fT6|) . Using Lemma 01 we get the estimate (fT7]). 



Q.E.D 



Remark. Using Sobolev embedding theorems, we get global existence of classical solutions 
provided that initial data are smooth. 



3 Biological background 

In order to study the general system Q - ©, we have to first specify realistic assumptions on 
the parameters of the model. To this end, we summarize the relevant biological processes in 
Section nill and we extract the mathematical assumptions in Sect ion [3.2 1 These assumptions will 
be later used to prove the global existence results in Section |11 

3.1 Bacterial chemotaxis 

As discussed before, the bacterial movement can be viewed as a biased random walk. Bacterial 
motility is commonly provided by flagella, which are long, spiral-shaped protein rods that stick 
out from the surface of the cell |2S] • The example of flagellated bacterium is the enteric bacterium 
E.coli which has 6-8 flagella. It has two modes of behaviour based on counterclockwise and 
clockwise flagellar rotation. When the flagella rotate counterclockwise (CCW), they all point in 
one direction and consequently the cell moves forward in a straight "run" . The speed of running 
is s = 10 — 20^m/sec. Clockwise (CW) rotation of the flagella causes the flagella to point in 
different directions, and the cell tumbles in place. Tumbling reorients the cell, so that it can 
move in new direction when running starts again. 

For E.coli, the duration of both runs and tumbles are exponentially distributed with means 
of 1 sec and 10~^ sec respectively if an extracellular chemical signal is not present 12]. Under the 
influence of an attractant, the cell increases its time in running in a favourable direction - see 
Figure ^ As the mean time for tumbling is ten times smaller than the mean time of running, we 
can often neglect the time spent tumbling and we can model the movement of the bacterium as 
a velocity jump process [201 ^1 1^ as we already did in Section ^ It means that the bacterium 
runs in some direction and at random instants of time it changes its direction with mean turning 
rate A(y). 

Since the bacteria move with more or less constant speed, the set V of all available velocities 
might be considered equal toV = sS'^~^ where is a unit sphere in R'^ and s is the speed of 

the bacterium. Let us note that set V = sS^~^ satisfies the general condition © (the presented 
theory works for any set V which satisfy ©). 
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(a) no attractant (b) increasing attractant (or decreasing repellent) 




Figure 1: (a) A typical bacterial trajectory when no attractant is present, (b) Under the influence 
of an attractant, the cell increases its time in running in a favourable direction. 



The kernel K{v,v' ,y) gives the probabihty of a change in velocity from v' to v, given that 
a reorientation occurs. The simples possibility is to assume that kernel is constant, i.e. 

i^(v,v',y) = p^. (32) 

This formula satisfies the normalization condition © . The underlying assumption behind (|32|) is 
that (during the tumble) bacterium simply choose a new direction randomly which is relatively a 
good approximation for the bacterial chemotaxis, although there is also some bias in the direction 
of the preceding run jH . More realistically, one can assume that the turning kernel is a function 
of the angle between new and old velocity, i.e. 

K{^r,^r',y) =k{e), where cos (0) = (33) 

Whatever the choice of K{v,v' , y) is, we may assume that it is bounded from above by a constant, 
i.e. 

i^(v,v',y)<C (34) 

where C is independent of v, v' and y. Next, we have to specify the choice of Q and the 
properties of the turning frequency A(y). 

Chemotaxis is the process by which a cell alters its movement in response to an extracellular 
chemical signal. From the microscopic (cell) point of view, bacterial chemotaxis consists of 
several steps. First, the cell detects the signal using its receptors. Then the signal information 
propagates through the signal transduction biochemical network described by ((SJ. The output 
of this network is a phosphorylated form of the protein CheY (denoted CheY-P) which alters the 
motor behaviour of the flagellar motors and consequently, the movement of the cell. CCW is the 
default state in the absence of CheY-P, which binds to motor proteins and increases CW rotation. 
Attractant binding to a receptor reduces the phosphorylation rate of CheY and thereby increases 
the time spent in running state which constitutes the fast response to a signal called excitation 
of signal transduction network. Another important aspect of signal transduction network is 
adaptation which means that the response (probability per unit time of CCW/CW rotation of 
flagella) returns to baseline levels on a time scale that is slow compared to excitation, provided 
that there is no further change in attractant concentration around the cell. 

A schematic of the signal transduction pathway is shown in Figure[21and it can be described as 
follows [271 1281 [T^ . Aspartate, the attractant most commonly-used in experiments (denoted S in 
FigurelSI), binds directly to the periplasmic domain of its receptor. Tar. The cytoplasmic domain 
of Tar forms a stable complex with the signaling proteins CheA and CheW (denoted A and W, 
respectively, in Figure |2)) , and the stability of this complex is not affected by ligand binding |14j . 
The signaling currency is in the form of phosphoryl groups (-P), made available to the CheY 
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(fast) excitation: attractant binding reduces 
autopliosphorylation rate of ClieA and 
it decreases concetration of ClieY-P 



(slow) adaptation: attractant binding increases 
metliylation state of receptor and it returns 
phospliorylation rates to baseline levels 




Figure 2: Excitation and adaptation in signal transduction pathway of E. coli (from 'AO] , with 
permission) . 

(denoted Y in Figure |2) and CheB (not shown in Figure \^ through autophosphorylation of 
CheA. Receptor complexes have two alternative signaling states. In the attractant-bound form, 
the receptor inhibits CheA autokinase activity; in the unliganded form, the receptor stimulates 
CheA activity. Consequently, the response of the signal transduction network to a step increase of 
the attractant concentration is as follows. First, the attractant binding to a receptor reduces the 
autophosphorylation rate of CheA. The level of phosphorylated CheA is thus lowered, causing less 
phosphate to be transferred to CheY, yielding a lowered level of CheY-P. As a result, tumbling 
is suppressed, and the cell's run length increases. This constitutes the excitation response of 
the system. Next slow methylation and demethylation reactions begin to influence the response. 
Ligand-bound receptors are more readily methylated than unliganded receptors, and the lowered 
level of CheA-P causes a decrease in the level of CheB-P, thereby reducing its demethylation 
activity. As a result, the equilibrium of the system shifts in the direction of the higher methylation 
states. The autophosphorylation rate of CheA is faster when the associated Tar-CheA-CheW 
complex is in a higher methylation state, and so there is finally a shift back toward the receptor 
states containing CheA-P. As a result, CheY-P returns to its prestimulus level, and thus so does 
the CW bias of the cell. This constitutes the adaptation response. These key steps, excitation 
via reduction in CheY-P, when a receptor is occupied, and adaptation via methylation of the 
receptors, have been already incorporated in the mathematical models of the bacterial signal 
transduction HZmiH. 

Since the turning rate of bacterium is altered by CheY we can write A(y) = A(yi) where 
y\ denotes the concentration of the phosphorylated form of CheY. Hence, the individual-based 
model for bacterial chemotaxis is fully specified by the equation Q which is integrated along the 
trajectory of each cell, and by the yi component of the solution together with A = \{yi). The 
essential aspects of the dynamics which must be captured by model (jH)) are (i) it must exhibit 
excitation, which here means a change in the turning frequency \{yi) in response to a stimulus, 
(ii) the bias must return to baseline levels (i.e., the response must adapt) on a time scale that 
is slow compared to excitation, and (iii) the signal transduction network should amplify signals 
appropriately [HH^ni- The mathematical assumptions on ^ and \{yi) are given in Section EUl 

3.2 Mathematical assumptions on the signal transduction network 

The mathematical model of the signal transduction network Q can be rewritten in the following 
form 



-^=F(C(t),y) 



where 




(35) 
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The vector function C(t) gives signal values which are seen by a cell along its trajectory. Time 
evolution of y in equation ()35p is controlled by the input time dependent vector C(t). Therefore, 
it is natural to describe the behaviour of F in terms of the input function C{t). 

The mathematical formulation of the adaptation property of the signal transduction network 
(jni) can be written in the following form. There exists a universal constant Iji such that for any 
constant signal along the trajectory Co, i.e. C(t) = Co = const and for any initial condition 
y(0) = yo, the solution of the system (jSJ satisfies 

limyi(t)=yi. (36) 

Formula (|36jl describes the perfect adaptation. From the application point of view, it is desirable 
that the signal transduction model satisfies (at least approximately) the adaptation property for 
a reasonably large set of signals. However, the existence theorems presented in Section |1] do not 
require perfect adaptation and we will prove the existence of solutions even for models which 
do not satisfy ()36() . It is worthwhile to note that the simplified model of excitation-adaptation 
dynamics Ullf) from Section |21 satisfied adaptation property In fact, yi(t) ^ as t ^ oo 

for any constant signal, i.e. yi adapts perfectly to any constant stimulus. Moreover, model 
describes the excitation-adaptation dynamics as discussed in Section ITT] provided that we choose 
te < ta- Here, the time constants te and ta are labeled in anticipation of using yi for the internal 
response, and 2/2 as the adaptation variable, and therefore we call te and ta the excitation and 
adaptation time constant, respectively jl2j . 

In order to model the random walk of the individual bacterium, we must have a good un- 
derstanding of the dependence of the (output) turning rate A(yi) on the (input) signal function 
C{t). If the input signal function is constant then the behaviour of \{yi) follows the adaptation 
property. On the other hand, time dependent input C{t) can introduce large variations in A(yi). 
The time derivative of C(i), i.e. the time derivative of the signal seen by a cell, is equal to 

dC dS dS 

To see what type of conditions on the turning rate A are reasonable, let us consider the time 
independent signal (attractant) with a maximum at the point Xm as it is schematically shown in 
one dimension in Figure 01 (panel in the middle) . We consider that bacteria move with the fixed 

density of extracellular density of 

bacteria signal bacteria 

(a) 



Figure 3: Schematic of behaviour of hypothetical cells which "perfectly avoid going in wrong 
directions" (panel on the left) and hypothetical cells which "perfectly follow good directions" (panel 
on the right). Details are explained in the text. 

speed either to the right or left and we discuss the following two simple cases of dependence of 
output A(yi) on input C(t). 

, , , , , f 1 for dC/dt > 0; „ ^ x , ^ f for dC/dt > 0; 

^^y^^ = I 00 for dC/dt < 0; ^^y'^ = I 1 for dC/dt < 0. 
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Let us note that cases (a) and (b) are considered as definitions of tlie input-output behaviour in 
two extreme cases (these definitions are not connected with any underlying differential equation 
in this example). 

First, suppose that a bacterium is at the position x < Xm- If we use input-output behaviour 
(a), then the cell goes to the right. It sometimes "turns" to the left but it instantly turns 
back. So, the cell spends all the time going to the right, and case (a) is an example of the 
individual-based model where cells perfectly avoid going in wrong directions. If we use input- 
output behaviour (b), then the right going cells never turn (for x < Xm)- Hence, case (b) is an 
example of the individual-based model where cells perfectly follow good directions. Both cases 
(a) and (b) describe the simple transport of bacteria for x < Xm- The difference of these models 
is when cells reach the maximum of the signal Xm- In case (a), cells instantly turn back. It means 
that the final positions of all bacteria are equal to Xm and a Dirac-like distribution is created in 
finite time (see Figure |31 panel on the left). In case (b), cells continue movement to the region 
X > Xm and the final distribution profile is smooth, as shown schematically in Figure |31 (panel 
on the right). 

The previous simple example shows that singularities might develop if the turning rate is too 
large (without a reasonable control by the signal change), as in case (a) where cells perfectly avoid 
going in wrong directions. This observation suggests for growth conditions on X{yi) from above 
which prevent formation of singularities. The necessary conditions on the turning frequency 
X{yi) is X{yi) > and our heuristic conclusions can be incorporated to the following growth 
estimate 

A(2/i)<C 1 + A(|C|) + 



dt 



(38) 



where A(-) G C (M) is a non-negative, nondecreasing continuous function. The verification of 
growth estimate depends on the particular form of F(-, •) and A(-). For example, if Q and 
A(-) satisfy 



|yi| <Ci 1 + 



dC 



dt 



), A(yi) <C2(l + |2/ir), u;a<l, 



(39) 



then H38() follows. There are several other conditions on F(-,-) and A(-) which also guarantee 
growth estimate (|38|) . Hence, we do not formulate our growth estimates in terms of F(-, •) and 
A(-), but we simply assume H38() directly in our existence theorems. Using formula ((T)), we can 
formulate the estimate (|38j) also in terms of the kernel T(v, v',y). 

Using estimate ()34() and definition ((7|, we can write the growth assumption on T in the 
following form 

T{v,v',y)<C\X{yi)\. (40) 

We also have to assume a growth assumption of Vy-F. In Theorem [21 we assume that there 
exists a non-negative, nondecreasing continuous function n(-) G C (M) satisfying 

|VyF(z,y)| <C(l + n(|z|)). (41) 

Notice that our simple model Hll|) satisfies H41|) . A different condition on Vy-F is studied also 
in Corollary E 
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3.3 Mathematical assumptions on the dynamics of the extracellular signals 



Various forms of R(S,n) can be considered. The simplest case from the mathematical point of 
view is when the extracellular signals are nutrients which are consumed by cells, i.e. 

R(S,n) = -KSn (42) 

where K is a diagonal nonnegative M x M matrix (with rate constants on the diagonal). One 
can also assume that the cells produce signals which are degraded at some rate, i.e. 

K{S,n) =n[ki,k2,...,kMf -KS (43) 

where ki, . . . , kM are rates of production of the different components of the signal and K is a 
diagonal nonnegative M x M matrix. If we allow the nondiagonal terms in matrix K, then the 
extracellular coupling of the signals (e.g. reactions between signals) is added to the model. One 
can also consider that some signals can be produced by cells and some signals can be degraded 
by cells, i.e. effectively combining (|42j) and (|43|) . Moreover, we can also assume that some signals 
can be attractants while other signals can be repellents etc. 

Depending on the model system, there are many possibilities to specify the dynamics of the 
extracellular signal. In what follows, we use (|43|) . However, it is possible to modify and prove 
the following existence theorems using different evolution equations for the extracellular signal 
too. The only requirement is that the evolution equation for the extracellular signal must satisfy 
suitable growth estimates similar to the estimates which are proven in Lemma IHl for H43|l. 



4 Global existence for the general signal transduction models 

In this section, we prove global existence results using the framework of Sections 13 . 21 and 13 .31 We 
will work in one-dimensional physical space, i.e. = 1 and we first assume the case of elliptic 
equations for the extracellular signals. Hence, system of equations © ~ reads as follows 

^ + V, • v/ + Vy • F(S(x), y)/ = T{v, v' , y) [f{v') - f{v)] dv' , (44) 

+ hn - kfS, = 0, i = l,...,M, (45) 

where di, ki and k^ are positive constants and n = n{x, t) is the macroscopic density of individuals 
at point x € M and time t given by (|1U|) . Now, we can formulate the existence theorem. 

Theorem 2 Let us assume ((23), and Assume that /o G n L°°(M x F x M™) and 

let initial condition Sq G [I^^'P(R)]*^ satisfies Then there exists a global solution of system 

^441 ^ d-^^D satisfying, for all t > 

/(•,-,-,t) G nL^(M X y X M™), (46) 

S{-,t) G [W^'P{R)]^, for alll<p< +oo, (47) 
and initial conditions /(•, •, •, 0) = /o(-, •, •) and S(-, 0) = So(-). 
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Remark. To avoid technicalities, we focus in Theorem [21 only on estimates of /. The 
results could be extended to W^'^ estimates under suitable growth assumptions on derivatives 
of T (v, v' , y) and F. 

In order to prove Theorem [21 we formulate some auxiliary lemmas. We start with the general- 
ization of the Gronwall inequality. 

Lemma 4 Let a{s) and b{s) be positive integrable functions on [0,t]. Let w (t) be positive and 
differentiable in t, and satisfy 

w' < a (t) wlnw + b (t) w. 



-1 exp 



(Jo <i{s)ds) 



Then 

w{t)< w (0) exp (^J^ b (s) e- '^(^)'^^ds 
Proof. See Lemma 4]. 

Q.E.D. 

The characteristics of the hyperbolic equation H44|) are given for = 1 as the solution of H18() . 
The back-in-time characteristics starting at (x, v, y, t) are given as 



X {s;x,v,y,t) = x — v {t — s) , 

Y{s;x,v,y,t) = y-[ F (S (X(r), r) , Y(r)) dr. 

J s 



(48) 
(49) 



The generalization of Lemma ^ is given as the following Lemma. 



Lemma 5 Derivation of the characteristics and with respect to the initial conditions 
gives 

dX , dY 



dx 



1 



and 



dy 



exp 



^ ^(S(X(r),r),Y(r))dr 



Moreover, 



, dY 

det -- — = exp 
dy 



^*VyF(S(X(r),r),Y(r))dr 



Proof. We differentiate (|49j) with respect to y to get 

dY 



dy + i ^(SWr),r),Y(r))-(r)dr. 

where Im is the m x m identity matrix. Let 

dV dY 

— (S(X(r),r),Y(r)) — (r)dr, 



(50) 
(51) 

(52) 



then we have 



dF dF 

G' (s) -G{s) — (S (X(r), r) , Y(r)) = — (S (X(r), r) , Y(r)) 
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Integrating the last equation, we obtain (|5U() . Since the determinant of the exponential of the 
matrix is the exponential of the trace of the matrix, we have 



, dY 

det — — = exp 
oy 



exp 



trace (^-^ ^ (S (X(t), r) , Y(r)) dr 
- rVyF(S(X(r),r),Y(T))dT 

J s 



Hence, we have proved ()51() . 



Q.E.D. 



Next, we present the growth estimates on the extracellular signal S and on its derivatives. The 
time and space derivatives of the signal vector S are controlled by logarithm of the L^-norm of 
the cell density. Note that the analogous result was also shown in |16| Lemma 4] for the parabolic 
equation for the extracellular signal. The difference between jl61 Lemma 4] and Lemma IHl is that 
we prove also the estimate on the time derivative as well as the estimate on the space derivative 
of the signal. 

Lemma 6 If n £ L~([0, oo) : (M) n (M)), then the solution S in (g^ satisfies 

\\SmL^<C\\n{t)\\L,=C\\n{0)\y, 



dS 


it) 


dx 


as 


it) 


'dt 



< 



< 



C[l + ||n(0)||^,{l + ln(||n(t)||^. + l)}] 



C[l + ||n(0)||ii{l + ln(||n(t)||^.+l)}]. 



(53) 
(54) 



where the constant C depends only on ki, , di and V . 

Proof. Let 1 < i < M. Taking the Fourier transform of H45|l in the x— variable, we obtain 

ki n (^, t) 



S^i^,t) 



di e + k^/d^ ■ 



Thus we have 



\Siit)\\ 



< 



S^it) 



, ki II ^ 

< — n 
Li di 



'h- fk- i-O' 



Next we estimate the x— derivative of the signal as follows 

< iSiit) 



\i\\nii.t)\ h , 



where 1\ 



\i\\nii.t)\ 



d^ and I2 



mnii,t)\ 

m\n{t)\\l2 C'^ + ki/d^ 



de- 



le 



First, we estimate the integral Ii. We obtain 

h < ||n(t)|Loc 



\n(t}\\^ac In 



+ 1 < ||n(t)||„ln 



\n{t)r 



L2 



+ 1 



We use Holder's inequality with p = q = 2 to estimate I2 as 



h < \\n 



< n 



L2 



L2 



\ 1/2 

r^de < V2. 



2 \ 1/2 

da 



By combining the estimates for Ii and we obtain H53|). In order to estimate the time derivative 
of the extracellular signal, we take the time derivative of ()45|) and apply the Fourier transform 
in the x— variable to get 

t) = h^(e t) - 

dt ^^''^ d,dt + 

By integrating (|44p over v and we get 

^TZ dj f f 

— = -— where j {x,t) = vf{x,v,y,t)dvdy. 
dt dx JJvxW^ 



Thus we have 



Then we have 



Notice that 



dt ^^'^> d^e + k^di' 



dSi 
dt 



< 



dt 



it) 



< 



ki 



lei ii(e)i 



< 



-1 < 




e+kVd, 

\v\ f {x, V, y, t) dxdt;dy < 

xyxR™ 



<C{V)\\n{t)\\^,=C{V)\\n{0)\\^,, 

(2 \ 1/2 
[ ( [[ \v\f{x,v,y,t)dvdy) dx] < 

<C{V)\\n[t)\\^, 

where we have used that V is compact. Using similar ideas as in the proof of estimate ()53|1. we 
prove 

Q.E.D. 
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Lemma 7 Let F satisfy l \41l - Then the characteristics (\4^ ^ satisfy for all < s < t, 



det — s 
dy 



< 



exp [Ct{l + n(C7||/o||iO}]- 



(55) 



Proof. Using Lemma El and (|4H) . we obtain 



det — s 

5y 



-1 



exp 



< exp 



< exp 



C [ l + n(|S(X(T),r)|)dT 

J s 

Using Lemma ini we deduce (|55p . 

Proof of Theorem 2. Using (|55)) and (|in)) . we obtain 

T{v,v',y)<c(l + Ai\C\) + 



^ VyF(S(X(r),r),Y(r))dT 

C7t|l + n( sup ||S(t 



Q.E.D. 



dC 



dt 

Integrating (|44|) along the characteristic (|48|) - (|49|) from to i and using (|56j) . we obtain 

/(2;,t;,y,t) </o(X(0),t;,Y(0)) + 



(56) 



A(|S|) + 



(X(r),r) X 



X y f {X{t),v',Y{t),t) dv'|dr+ 

+ r I Vy • F {S{X (r)), Y (r)) \f {X (r) , t;, Y (r) , r) dr, 
Jo 

where we used that V is compact. By virtue of assumption H41|) . |Vy • F| is bounded by 
C7(l + n(|S(X(r),r)|)). Thus we have 



f{x,v,y,t)<foiX{0),v,Y{0)) + 



(57) 



A(|S|) + 



dS 



dt 



+ 



dS 



dx 



(X(r),r) X 



X J^f{XiT),v',YiT),T) di;'|dr+ 

+C7(y) /* (1 + n (IS (X(t), r) D) / (X (r) , V, Y (r) , r) dr. 
Jo 

Using Lemma |71 we obtain for t > 0, 

[ [ fP{X{T),v',Y{T),T)dv'dxdvdy = 

jRxVxRM Jy 

= \V\j f^X (r) , v', Y (r) , r) (det ' (det 



dX 

dx 



dv'dXdY < 
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< \V\ e^* j fP {X (r) , v', Y (r) , r) dv'dXdY. 
We take the p-th power of (jFTj) and integrate over x, v, and y to get for t > 0, 



ll/(t)llLP<ll/ollL.+C(^) 



dS 



+ 



t 



dS 

dx 



;i + A(||S(r)|L^) + n(||S(r)||^^)+ (58) 



X 11/ W IlLpdr. 



Using Lemma IHl we get for all < t < T and with p = 2 



\\f{t)h. < Wfoh^ + Ce^l [1 + ||n (0)11^1 {1 + ln(||n(r)||^. + 1)}] x ||/(r)||^,}dr < 

Jo 

< Wfoh. + Ce^'f [1 + 11/ (O)ILi {1 + In (||/(r)||^. + 1)}] x ||/(r)||^,}dr. 
Jo 

By applying the Gronwall Lemma to 1)59^. we obtain for t >0, 

WfmL^ <C {h,k^,d,,A,U,V,\\fo\\Li,\\foh2) < 



(59) 



oo, 



— (t) 

dx 



<C {h,k^,di,A,U,V,\\fo\\LiA\fo\\L2) < oo, 



(60) 



where we used Lemma IHl to get estimate (jHOl- We now apply (|6fl]) to (|58|) and we get (for t >0 
and for all 1 < p < oo) 



||/(t)|lLP <c(fc„fcO,d„ A, n,y,||/o||^,,||/o||^„ 11/011^2) < 



00, 



(61) 



i.e. we have obtained Using the elliptic equation H45() . the second derivative of the extra- 

cellular signal can be expressed as 



dx"^ dj di 



(62) 



Using l|6T1) and the elliptic theory, we deduce (|17]) . Thus we complete the proof of Theorem |2j 

Q.E.D. 

We conclude this section with two corollaries. They provide other conditions for the global 
existence of solutions. The proofs are omitted because they are similar to proofs of Theorem ^ 
and Theorem [21 

Corollary 1 Assume ([23) and ^4(^ - Suppose there exists a non-negative, nondecreasing contin- 
uous function 11 (•) G C (M) and 7 > with ujj < 1 satisfying 

VyF<0 and VyF (z, y) < C (1 + U (|z|) + |y|^) (63) 

Assume that /o G n L~(M xVx M™) and let the initial condition Sq € [Ty2'P(M)]^'^ satisfy 
0451 ). Then there exists a global solution of the system l \44l ~ satisfying, for all t > 

/(•, •, ■,t) eL^n L°°(M xV X M™), 

S(-, t) G [W^'P{R)f\ for alll<p< +00 

and initial conditions /(•, •, -,0) = /o(-, •, •) and S(-, 0) = So(-)- 
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Corollary 2 Assume 

A(yi) < C, Tiv, v', y) < C{1 + I A(yi)l). (64) 

We further assume that F satisfies either (jTTJ or ^U^. Assume that /o G n L°°(]R xV x M™-) 
and So G [VF^'°°(R)]*^ wzi/i compact support. Then there exists a global solution of system of 
equations H44\) \14\ satisfying 

f{-,-,-,t) G nL°°(M xV X M™), 

S(-,t) G [W^i'°°(M)]*^ 
anrf initial conditions /(•, •, •, 0) = /o(-, •, •) and S(-, 0) = So(-)- 

5 Discussion 

The simplified model of the bacterial signal transduction was studied in |12| where equation 
© was given as ((TT|) . Using model (|TT|) for the steady extracellular signal, one can derive the 
closed macroscopic ( Keller- Segel, chemotaxis) equation for some parameter regimes. See 
in ID and ^Hj in 2D/3D. Hence, the transport equation framework can be used to study the 
macroscopic behaviour in terms of microscopic parameters for the steady extracellular signals 
and simplified models of the signal transduction. 

Here, we focused on more complex models where we coupled the complex transport equation 
© with the parabolic or elliptic equation for the signal ©• The starting point of the analysis of 
such complex models is the existence theory. In this paper, we provided several sets of sufficient 
conditions for the global existence of solutions of system Q ~ ©• There are many open questions 
remaining, e.g. the existence theory in A^-dimensional physical space. It is also not clear whether 
one can derive the closed evolution equation for the density of cells n{x,t) as we did for the simple 
case of noninteracting particles jl2l I13j . If we are not able to derive the macroscopic equations 
then suitable computational approaches have to be used to study the macroscopic behaviour of 
bacteria [llj . 

There are several related results on kinetic models of the cellular movement. They often do 
not take the intracellular dynamics into account. Kinetic models were derived in ^ |^ using 
stochastic models of the movement of cells like bacteria or leukocytes. Reference |21j addresses 
the formal diffusion limit of kinetic models to the classical Keller-Segel model. The discussion 
on issues of aggregation, blow-up, and collapse for certain class of random walks can be found 
in |2n]- A Boltzmann-type kinetic model for chemotaxis without the internal dynamics coupled 
with an elliptic equation for the extracellular signal is studied in [7j where global existence and 
rigorous diffusion limit to the Keller-Segel model were proven. In jl5|ll7j. a more general kinetic 
model was treated in two and three dimensions. A one-dimensional hyperbolic model was studied 
in The papers |E1 El EI took into account the effect of the gradient and the temporal 

derivative of the chemical signal and showed the global existence of smooth solutions with smooth 
initial data as well as the rigorous diffusive limit to the classical Keller-Segel model. However, 
all the rigorous global existence results so far have not included the temporal derivative of the 
signal in the growth condition of the turning frequency as we did in this paper. See also |24j for 
more related works. 
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